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We examine the relation between two approaches to the quantum relativis- 
tic two-body problem: (1) the Breit equation, and (2) the two-body Dirac 
equations derived from constraint dynamics. In applications to quantum elec- 
trodynamics, the former equation becomes pathological if certain interaction 
terms are not treated as perturbations. The difficulty comes from singularities 
which appear at finite separations r in the reduced set of coupled equations 
for attractive potentials even when the potentials themselves are not singular 
there. They are known to give rise to unphysical bound states and resonances. 
In contrast, the two-body Dirac equations of constraint dynamics do not have 
these pathologies in many nonperturbative treatments. To understand these 
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marked differences we first express these contraint equations, which have an 
"external potential" form similar to coupled one-body Dirac equations, in a 
hyperbolic form. These coupled equations are then re-cast into two equiv- 
alent equations: (1) a covariant Breit-like equation with potentials that are 
exponential functions of certain "generator" functions, and (2) a covariant or- 
thogonality constraint on the relative momentum. This reduction enables us 
to show in a transparent way that finite-r singularities do not appear as long 
as the the exponential structure is not tampered with and the exponential 
generators of the interaction are themselves nonsingular for finite r. These 
Dirac or Breit equations, free of the structural singularities which plague the 
usual Breit equation, can then be used safely under all circumstances, encom- 
passing numerous applications in the fields of particle, nuclear, and atomic 
physics which involve highly relativistic and strong binding configurations. 
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I. INTRODUCTION 



In contrast to the accepted description 

[7 -p + m + V]ip = (1.1) 

of the relativistic quantum mechanics of a single spin-one-half particle moving in an exter- 
nal potential V given by Dirac, a number of different approaches have been used for two 
interacting spin-one-half particles. A traditional one based on the Breit equation [1], also 
known as the Kemmer- Fermi- Yang equation [2], 

[«i • Pi + Ami + a 2 ■ p 2 + ftms + V{r 12 )}ip = Eip, (1.2) 

contains a sum of single-particle Hamiltonians and an interaction term between them. (E 
is the total energy in an arbitrary frame.) Although the Breit equation is not manifestly 
covariant, it has provided good perturbative descriptions of the positronium and muonium 
energy levels. However, it is well known that those parts of the Breit interaction 

-a I 1 

V(r) = (1 - -ai • ct 2 - -cki • ra 2 ■ r), (1.3) 

r 2 2 

in the Breit equation beyond the Coulomb term should not be treated nonperturbatively, 
but must be handled only perturbatively. In other words, a consistent treatment of the Breit 
equation in powers of a will generate unwanted terms not present in quantum electrody- 
namics [Refs. 3,4]. 

A more serious difficulty of the Breit equation is that when the interactions are treated 
nonperturbatively, structural pole singularities could appear at finite r even when the in- 
teractions themselves are singularity-free there [5,8]. Two of us have found recently that 
these pole singularities occur under certain conditions depending on well-defined algebraic 
relations among the different potentials that could appear ||[7] . They lead to unphysical 
states or unphysical resonances and therefore must be strictly avoided 

The primary purpose of this paper is to show how these pole singularities can be avoided 
from the beginning so that the Breit equation can be used without difficulty in diverse 
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applications in particle, nuclear, and atomic physics involving highly relativistic motions 
and strong binding potentials. This is accomplished by relating this older approach to one 
that has been developed much more recently. 

Dirac's constraint Hamiltonian dynamics [9] provides a framework for an approach pro- 
posed by Crater and Van Alstine [10,11] that differs notably from that of the Breit equation. 
It gives two-body coupled Dirac equations, each of which has an "external potential" form 
similar to the one-body Dirac equation in that for four-vector and scalar interactions one 
has 

= 751 [7i • (Pi - M) + = 0, (1.4a) 

S 2 ifj = 7 52 [72 • (P2 - A 2 ) + m 2 + S 2 }4> = 0, (1.4b) 

(ibi, i — 1,2 are the 75 matrices for the constituent particles). Unlike the Breit approach, 
these equations are manifestly covariant and have interactions introduced by minimal sub- 
stitution. 

They have common solutions if the operators Si,S 2 commute. This situation, called 
strong compatibility [9], can be achieved by the proper choice of the operators Ai and 
Si [10,11]. The commutator cannot, however, be made to vanish for more general types of 
interactions such as pseudoscalar or pseudovector. Nevertheless, under certain circumstances 
the commutator can be reduced to a combination of Si and S 2 themselves. The equations 
are then said to be weakly compatible [9], because this will also ensure that solutions of 
Slip = in the more general cases could be solutions of S 2 ip = as well. These compatibility 
properties are important because they guarantee the existence of common solutions to (i.e. 
the consistency of) the constraint equations before they are actually solved. 

Although the constraint two-body (CTB) Dirac equations have been used far less fre- 
quently than the Breit equation, they have important advantages over the latter. In de- 
scribing electromagnetic bound states [11,12] they yield nonperturbative and perturbative 
results in agreement with each other. That is, the exact (or numerical) solution produces 
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a spectrum that agrees through order a 4 with that given by perturbative treatment of the 
Darwin, spin-orbit, spin-spin, and tensor terms obtained from the Pauli reduction. In par- 
ticular, total cm. energies w for the e + e~ system in the 1 Jj states is found to satisfy a 
Sommerfeld formula [11,12] 



w = m 



2+2/1+ 



a 2 



\ I V n + .J(J + \) 2 -a 2 - J- ±) 2 

ma 2 ma A lima' 1 fi . . „ s 

= 2m - - + — + 0(a 6 ). 1.5 

An 2 2n 3 (2J+l) 64 n 4 1 ; v ; 

They agree through order a 4 with those of the perturbative solution of the same equation, 
and also with those of standard approaches to QED. A recent paper has numerically extended 
this agreement at least to the n — 1,2, 3 levels for all allowable J and unequal masses fl2"fl . 

In this paper, we are concerned with another advantage of the CTB Dirac equations, 
namely that no unphysical states and resonances of the type discussed in [6,7] have ever 
appeared in past applications. We are able to show here that this is in fact true for a general 
interaction and that this is a consequence of the exponential structure of the interactions 
appearing in them. This result is obtained by first reducing the CTB Dirac equations 
to a Breit equation and an equation describing an orthogonality constraint. The equivalent 
Breit equation is then shown to be singularity-free, provided that the exponential interaction 
structure is not destroyed by inadvertent approximations and that the operators appearing 
in the exponent are themselves free of finite r singularities. 

The exponential structure that tames the unphysical singularities turns out to be a 
consequence of a relativistic "third law" describing the full recoil effects between the two 
interacting particles. We carefully trace, in the formulation of the CTB Dirac equations, 
how this structure arises from the need to make these equations at least weakly compatible. 
Compared to the laissez-aller approach of the Breit equation for which any interaction 
seems possible, the restriction of the interaction structure needed in the constraint approach 
represents a conceptual advance in the problem. For this reason, we take pain to elucidate 
its conceptual foundation as we present elements of the CTB Dirac equations needed for our 
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demonstration that they are singularity-free. 

This paper is organized as follows. We start in Sec. II with a brief overview of the 
derivation of the constraint two-body Dirac equations for scalar interactions both to define 
the notation used and to describe the main concepts involved. One of the most important 
properties is the compatibility of the two constraints. We remind the reader in Sec. Ill 
how to introduce general interactions into them that preserve this property. In Sec. IV we 
derive a covariant version of the Breit equation from the constraint equations, introducing 
the concept of exponential generators for a wide range of covariant interactions. In Sec. V 
we clarify the structure of this covariant Breit equation by decomposing it into a matrix form 
involving singlet and triplet components of the matrix wave function followed by reduction 
to radial form by the use of a vector spherical harmonic decomposition. This reveals clearly 
how the constraint approach avoids the structural pole singularities that have plagued the 
original Breit equation since its introduction over 65 years ago. In Sec. VI we show by 
contrast, how the pole singularities arises for each of the nonzero generators if one uses the 
original Breit interaction. Section VII contains brief concluding remarks. 

II. TWO-BODY DIRAC EQUATIONS OF CONSTRAINT DYNAMICS 

Following Todorov [13], we shall use the following dynamical and kinematical variables 
for the constraint description of the relativistic two-body problem: 

i. ) relative position, x\ — x 2 

ii. ) relative momentum, p = (e2Pi — e\P2)/w , 



and 



iii. ) total cm. energy, w = \/—P 2 , 

iv. ) total momentum, P = p\ + p 2 , 

v. ) constituent on-shell cm. energies, 



w 2 + ml — m\ 



£2 = 



w 2 + m\ — m\ 



(2.1) 



ei = 



2w 



2w 



In terms of these variables, we have 
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p 1 = e 1 P+p, p 2 = e 2 P-p, (2.2) 

where P = P/w. 

We start from the (compatible) Dirac equations for two free particles 

<W> = (0i-Pi + ™Ai)V> = O, (2.3a) 

S 2 o*P = (02 ■ P2 + m 2 52 )^ = , (2.3b) 

where ip is just the product of the two single-particle Dirac wave functions. These equations 
can be written as 

S 10 ^ = (0i • P + tidi ■ P + mAOV' = , (2.4a) 

S 20 ^ = (-02 ■ V + £202 ■ P + m 2 52 )^ = 0, (2.4b) 

when expressed in terms of the Todorov variables. The "theta" matrices 

0? = ^ = 0,1,2,3,2 = 1,2 (2.5) 

05i = i \^ r r5i (2.6) 
satisfy the fundamental anticommutation relations 

[or]+ = -<r, (2.7) 

[0 5 *,0f]+ = O, (2.8) 

[0 5i ,0 5i ] + = -l. (2.9) 
[Projected "theta" matrices then satisfy 

[0 i -P,0 i -P]+ = l, (2.10) 
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[0,-P,0£l]+ = O, (2.11) 



where 9f ± = Qi V {j] llv + P^P")]. They are modified Dirac matrices 0] which ensure that the 
Dirac operators <Si and S 2 o ar e the square root operators of the corresponding mass-shell 
operators —\{p\ + rn\) and — \{p\ + m|). 

Using the Todorov variables and the above brackets, the difference 

(S 2 l0 - S 2 2Q )i, = = ^(pl + mi-pj- mj)^ (2.12) 

leads to an equation 

p .pip = -[ w ( ei - e 2 ) - (m? - m^)]V = 0. (2.13) 

The physical significance of the orthogonality of these two momenta is to put a constraint 
on the relative momentum ( eliminating the relative energy in the cm. frame). 

We will also use covariant (cm. projected) versions of the Dirac a and (3 matrices here 
defined by 

p i = - 7i .p = 26 5i 6i ■ P, (2.14) 
af = 2^-P, (2.15) 

and 

erf = TBiQ £ = 2V2i9 5i 6i ■ P6 ±i , i = l,2. (2.16) 

These covariant Dirac matrices take on the simple form af = (0, CKj) and erf = (0, cr,) in the 
center-of-mass system for which P = (1,0). 

If we now introduce scalar interactions between these particles by naively making the 
minimal substitutions 

rm -> Mi = rrii + Si, i = 1,2 (2.17) 
(as done in the one-body equation), the resulting Dirac equations 



SiV = (0i'P + e\Oi ■ P + = 0, 



(2.18a) 



S 2 ?P = {-6 2 ■ V + £202 ■ P + M 2 9 52 )^ = (2.18b) 
will not be compatible because 

ft, s 2 ]_V = [0i • p, M 2 e 52 ] + [m^ 51 , -e 2 - P }- = -i(dM 1 ■ e 1 e 52 + om 2 ■ e 2 e 51 )^ ± 0, (2.19) 

where d is the four-gradient. 

In the earlier work ||10"|,15|, compatibility is reinstated by generalizing the naive S\ and S 2 



operators with the help of supersymmetry arguments. The procedure contains four major 
steps: 

a) Find the supersymmetries of the pseudoclassical limit of an ordinary free one-body 
Dirac equation. 

b) Introduce interactions of a single Dirac particle with external potentials that preserve 
these supersymmetries. For scalar interactions, this requires the coordinate replacement 

9^9 

x^x^ = x^ + i . (2.20) 

(Note that the Grassmann variables satisfy 9 2 = 0. As a result this self referent or recursive 
relation has a terminating Taylor expansion). 

c) Maintain the one-body supersymmetries for each spinning particle through the re- 
placement 

(xi - x 2 ) -> (xi - x 2 ) (2.21) 

in the relativistic potentials S^. 

These steps lead to the pseudoclassical constraints (the weak equality sign « means these 
equations are constraints imposed on the dynamical variables) 

S x = (0! • v + eiOi ■ P + M x 9 hl - idM 2 /M x ■ 9 2 9 52 9 51 ) « 0, (2.22a) 
S 2 = (-0 2 ■ p + e 2 9 2 ■ P + M 2 9 52 + tdM 1 /M 2 ■ 9,9 52 9 51 ) « 0. (2.22b) 
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They are strongly compatible under the following two conditions: 

i.) the mass potentials are related through a relativistic "third law" 

d(Ml - Ml) = 0, (2.23) 

and ii.) they depend on the separation variable only through the space-like projection 
perpendicular to the total momentum 

Mi = Mi(x ± ), (2.24) 

where 

a*=(rr + P' i P v )(xi-x 2 ) v . (2.25) 
Integration of the "third law" condition yields 

Mj - Ml = m\-m\ (2.26) 

with the hyperbolic solution 

Mi = mi coshL + m 2 sinhL, M 2 = m 2 coshL + m\ sinhL, (2.27) 

given in terms of a single invariant function L = L(x±). 

The x± dependence of the potential and the relativistic "third law" lie at the heart 
of two-body constraint dynamics. Without these conditions the constraints would not be 
compatible. While the physical importance of the x± dependence lies in its exclusion of the 
relative time in the cm. frame, the "third law" condition relates the mutual interactions 
between the particles to the effective potentials each particle feels in the presence of the 
other in a consistent way. It is useful to show its implications in the simpler case of spinless 
particles. The two generalized mass shell constraints that are the counterparts of Eq. (2.22) 
for scalar interactions are 

Wi=p? + M?«0, 2 = 1,2. (2.28) 

The compatibility condition for these two constraints involves the classical Poisson bracket 
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[Hi,H 2 ] = [plMl] + [MIpI] + [Ml Ml] « 0. (2.29) 

One can see that this is satisfied provided that the "third law" condition Eq. (2.23) and 
condition (2.24) are satisfied. (Although the "third law" solution Eq. (2.26) combined with 
Eq. (2.24) is the simplest solution, it is not unique [|T5| .) 
For scalar interactions parametrized by 

Mi = m; + Si, z = l,2, (2.30) 

the "third law" condition becomes 

mi Si = m 2 S 2 , (2-31) 

in the nonrelativistic limit (|Sj| << m.j). This result can also be obtained from Eq. (2.27) 
by keeping only terms linear in L. The two constraints (2.28) can now be written as 

p 2 + M? « p 2 + 2miSi + S 2 + m 2 -e 2 = 0, (2.32) 

where we have used the fact that Hi — Ti 2 = P • p ~ remains unchanged upon the 
introduction of scalar interaction in Eq. (2.24). Hence the total cm. energy w = e± + e 2 
takes on a familiar form in the nonrelativistic limit 

2 

w = mi +m 2 + ^- + S + 0(S 2 ), (2.33) 
2/i 

where 

S = (nil + m 2 )S 1 /m 2 = (mi + m 2 )S 2 /mx. (2.34) 

d) The final step is to canonically quantize the classical dynamical system defined by S\ 
and c>2 by replacing the Grassmann variables 6^,6^, i — 1, 2 with two mutually commuting 
sets of theta matrices, and the position and coordinate variables by operators satisfying the 
fundamental commutation relation 

{x»,p u } -> [x^,p u ] =irf v . (2.35) 



11 



The compatible pseudoclassical spin constraints S\ and S 2 then become commuting quan- 
tum operators 

[«$!,&]= 0. (2.36) 
The resulting CTB Dirac equations for scalar interactions 

Slip = (9 1 -p + e 1 9 1 ■ P + Mi0 5 i - idL ■ 9 2 9 52 9 51 )ij = 0, (2.37a) 

S 2 ip = (~9 2 ■ V + e 2 9 2 • P + M 2 9 52 + idL ■ 9^9^ = 0, (2.37b) 

where 

dM 1 dM 2 

dL = -T7 1 = "TTi 2.38 
M 2 Mi 

are then said to be strongly compatible. This strong compatibility has been achieved by a 
supersymmetry which produces the extra spin-dependent recoil terms involving dL. These 
extra terms vanish, however, when one of the particles becomes infinitely massive (as seen 
by the parametrization Mj = m 8 + Si of the scalar potential) so that we recover the expected 
one-body Dirac equation in an external scalar potential. 
Note that the Dirac constraint operators satisfy [10] 

Si - S 2 2 = ~(pl +ml- P l-m 2 2 ) = -P-pnO. (2.39) 

Thus the relative momentum remains orthogonal to the total momentum after the intro- 
duction of the interaction. This also implies that the constituent on-shell cm. energies q 
are weakly equal to their off mass shell values (— P ■ pi q). Notice further that since 
[P-p, M(x±)] ~ P-x± = 0, this constraint does not violate the requirement of compatibility 
given in Eqs. (2.24-25). 

III. A GENERAL INTERACTION FOR TWO-BODY DIRAC EQUATIONS 

The previous work [16] has shown how the compatibility problem can be solved without 
having to invent new supersymmetries if the scalar potential is replaced by vector, pseu- 
doscalar, pseudovector, or tensor potentials. That work also relates the supersymmetric or 
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"external potential" approach to the alternative treatment of the two body Dirac equations 
of constraint dynamics presented by H. Sazdjian [17]. 

The "external potential" form Eqs. (2.37) of the CTB Dirac equations for scalar inter- 
action can be rewritten in the hyperbolic form [16] 

Slip = (coshA Si + sinh A S 2 )ip = 0, (3.1a) 

S 2 ip = (cosh A S 2 + sinh A Si)^ = 0, (3.1b) 
where A generates the scalar potential terms in (2.37) provided that 

A = -6 51 6 52 L(x ± ). (3.2) 
The operators Si and S 2 are auxiliary constraints of the form 

SiV> = (Sio cosh A + S 20 sinh A )ip = 0, (3.3a) 

S 2 ip = {S 20 cosh A + 5 W sinh A = 0. (3.3b) 

To verify that the "external potential" forms Eq. (2.37) result from using Eqs. (3.3) in 
Eqs. (3.1), one simply commutes the free Dirac operator S i0 to the right onto the wave 
function using Eqs. (2.7-2.11), (2.38) and hyperbolic identities [16]. With this construction, 
the interaction enters only through an invariant matrix function A with all other spin- 
dependence a consequence of the factors contained in the kinetic free Dirac operators S\o 
and iS 2 o- 

Even though the form of the contraints Eqs. (3.1) and (3.3) were motivated by exam- 
ining world scalar interactions, let us propose them for arbitrary A and determine their 
compatibility requirements. We do this for arbitrary interactions by generalizing arguments 
given in Refs. [16-17]. First consider the conditions for the compatibility of Eqs. (3.3a-b). 
Multiplying Eq. (3.3a) by <Si and Eq. (3.3b) by S 20 and subtracting we obtain 

P-p (cosh A) if> = 0. (3.4a) 
13 



Multiplying Eq. (3.3b) by <Sio and Eq. (3.3a) by S20 an d subtracting we obtain 

P-p (sinh A) i) = 0. (3.4b) 

We have used Eq. (2.2) and e x — e 2 = (m 2 —m\)/w to simplify these equations. Multiplying 
Eq. (3.4a) by sinhA, Eq. (3.4b) by coshA, bringing the operator P ■ p to the right and 
subtracting we find the condition 

[P • p, A]-0 = 0. (3.5) 

Multiplying Eq. (3.4a) by coshA, Eq. (3.4b) by sinhA, bringing the operator P ■ p to the 
right and subtracting we find the further condition 

P-p?p = 0. (3.6) 

Notice that this latter condition was previously associated with the "third law" condition 
when derived from the "external potential" forms of the constraints (see Eq. (2.39)). Here 
the "third law" condition is built into the constraint by having the same generator A for Eqs. 
(3.3a) and (3.3b). Thus the two tentative constraints Eqs. (3.3a) and (3.3b) taken together 
imply that for arbitrary A the orthogonality condition P ■ p pa has to be satisfied when 
acting on ip. However, in order to verify that there are no additional conditions beyond Eqs. 
(3.5) and (3.6) we must check for mathematical consistency by examining the compatibility 
condition. We compute the commutator [Si, S 2 ] by rearranging its eight terms and find that 

[Si,S 2 ] = [Si ,coshA]S 2 - [S 20 ,coshA]Si + [<S 20 , sinhA] S 2 - [«Si , sinhA]Si 

+coshA(S 2 - S 2 2 ) sinhA - sinhA (S 2 - S 2 ) coshA (3.7) 

does not in general vanish, unlike Eqs. (2.36) and (2.37). By using Eqs. (2.12) and bringing 
the operator P -p to the right, and using the conditions given in Eqs. (3.5) and (3.6) we can 
reduce [Si, S 2 }^> to only terms involving Slip and S 2 ^. Since Eqs. (3.5) and (3.6) follow from 
combining the constraints S^ = 0, no further conditions for mathematical consistency need 
be imposed on the constraints or the wave function. Eq. (3.6) is the quantum counterpart of 
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Eq. (2.39) but for arbitrary interactions. Eq. (3.5) is also satisfied for arbitrary A provided 
only that the generator A satisfies 

A = A(x±) (3.8) 

generalizing Eq. (2.24). 

The weak compatibility of the "external potential" form of the constraints Eq. (3.1) for 
general A 

[«Si,S 2 W> = (3.9) 

can be seen by examining the four commutators in [iSi,iS 2 ]. The commutator 

[coshA Si, coshA S 2 ] = coshA(coshA[S 1 , S 2 ] + [Si, coshA]S 2 + [coshA, S 2 ]Si) (3.10) 

and is weakly zero since S^ = and [Si, S 2 ]^ = 0. Likewise, we have 

[sinhAS 2 , sinhASi]^ = 0. (3.11) 

The remaining two brackets are 

- [coshA Si,sinhA Si] - [sinhA S 2 , coshA S 2 ] 
= [Si,sinh(A)]Si -sinh(A)[cosh(A),Si]Si + (1 -> 2). (3.12) 

Since they contain the constraints on the right we obtain Eq. (3.9) after combining with 
Eqs. (3.10) and (3.11). 

One final feature should be mentioned. Eqs. (3.1) and (3.3) are also applicable for a sum 
of the four "polar" interactions 

Ai = Al + Aj + Ag + A r , (3. 13) 

where 

A L = -L9 51 9 52 scalar, (3.14) 
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Aj = JP ■ 9 X P ■ 9 2 time - like vector, 



(3.15) 



Ag = Q9i± ■ 9 2 ± space — like vector, (3.16) 

A r = 4F9 1± ■ 9 2± 9 51 9 52 P ■ 9 X P ■ 9 2 tensor (polar). (3.17) 

For the sum 

A 2 = A c + A H + A/ + A y (3.18) 

of their axial counterparts 

Ac = — C/2, pseudoscalar, (3.19) 

A H = -2HP ■ 9 X P ■ 9 2 9 51 9 52 time - like pseudovector, (3.20) 

A/ = 2I9i± ■ 9 2 ±9 5 i9 52 space — like pseudovector, (3-21) 

A y = -2Y9 1± ■ 9 2± P ■ 9 X P ■ 9 2 tensor (axial), (3.22) 

the sinh A terms in Eq. (3.1) should carry negative signs instead [16]. In contrast, Eq. (3.3) 
as written remains valid as is for A 2 . For systems with both polar and axial interactions 
[16], one uses A x - A 2 in (3.1) and Ai + A 2 in (3.3). The terms L, J, Q, J 7 , C, H, I, Y are 
arbitrary invariant functions of x±. 

Eq. (3.3) is more convenient to use for the construction of Breit-like equations from the 
constraint formalism for general interactions. Eq. (3.1) is more convenient if one aims to 
obtain a set of Dirac equations in an "external potential" form similar to that exhibited in 
the one body Dirac equation with the transformation properties one would expect for such 
potentials. We have already seen this for the scalar case in which the scalar interaction 
"generator" L in (3.14) in the hyperbolic form leads to a modification of the mass term. In 
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the case of combined scalar, time- and space-like vector and pseudoscalar interactions, we 
use the hyperbolic parametrization 

Mi = mi coshL + m 2 sinhL , (3.23a) 
M 2 = m 2 coshL + mi sinhL , (3.23b) 



Ei = €i cosh J +e 2 sinhJ, (3.24a) 
E 2 = e 2 cosh J + e 1 sinhJ, (3.24b) 

G = e G , (3.25) 

where L, J, and Q generate scalar, time-like vector and space-like vector interactions respec- 
tively, while C generates pseudoscalar interactions. The resultant "external potential" form 
for 

A = Aj + A L + A g + A c (3.26) 

is 

Siij = (G9i ■ p + EiOi ■ P + Mi0 5 i 

+iG{6 2 ■ dQ9 1± ■ 9 2± + 9 2 ■ dJdi ■P9 2 -P-9 2 - dL9 51 9 52 + 9 2 ■ dC/2))^ = (3.27a) 

S 2 ^j = ( - G9 2 ■ p + E 2 9 2 ■ P + M 2 9 52 

-iG(9 1 ■ 0G9i ± ■ 9 2± + 9i ■ dJ9 x ■ P9 2 P - 9 X ■ 0L9 5 i9 52 + 9 X ■ dC/2))$ = 0. (3.27b) 

The scalar generator produces the mass or scalar potential M; terms, the time-like vector 
generator produces the energy or time-like potential Ei terms, the space-like vector generator 
produces the transverse or space-like momentum G terms, while the pseudoscalar generator 
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produces only spin-dependent terms. Note that the vector and scalar interactions also 
have additional spin-dependent recoil terms essential for compatibility. The above features 
are just what one would expect for interactions transforming in this way and are direct 
consequences of the hyperbolic parametrization of the contraints. Other parameterizations 
may not have this important property. 

IV. REDUCTION TO A BREIT EQUATION 

We now derive a Breit equation from the CTB Dirac equations (3.3). Consider the 
combination 2(6*i ■ PSi + 6 2 ■ PS 2 ). The terms proportional to w = e\ + e 2 are 

u? [cosh A + 2(0i • PQ 2 ■ P) sinh A] = w expiV) (4.1) 

where 

V = 2(6 l -P 9 2 -P)A. (4.2) 
The other terms can also be written in terms of V, using the relations 

cosh £> = cosh A, sinh V = 28 1 ■ P9 2 ■ P sinh A, (4.3) 

(26 2 ■P6 1 -p- 26 l ■ P9 2 ■ p) sinh A = -(26 1 ■Pd 1 -p- 26 2 ■ P6 2 ■ p) sinh£>, (4.4) 

and 

(26 2 ■ P0 5 i7Tii + 2Q X ■ P6 52 m 2 ) sinh A = -(29 l ■ P6 51 m 1 + 26 2 ■ P6 52 m 2 ) sinh V. (4.5) 

Then the combination 2(8 1 ■ PSi + 8 2 ■ PS 2 ) of the two hyperbolic equations in Eq. (3.3) 
takes the simple form 

w exp{V)^ = (H w + H 20 )exp(-V)ij (4.6) 

where 
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if io = -201 ■P9 1 -p- 20i • P0 5 imi = ai -p± + Ami + P ■ p 



(4.7) 



H 20 = 20 2 • P0 2 • p - 20 2 • P6 52 m 2 = -a 2 ■ p ± + (3 2 m 2 -P-p (4.8) 

(with the definition p± = p + Pp ■ P) are covariant free Dirac Hamiltonians involving the 
covariant a and f3 matrices given previously in Eqs. (2.14) and (2.15). If we take 

= exp{-V)i) (4.9) 

we obtain finally the manifestly covariant Breit-like equation 

w exp(2V)^ = (#io + H 20 )^. (4.10) 

The interactions appear in the Breit equation in the exponentiated form exp(2V), where 
TD contains all eight interactions shown in Ai and A 2 . We can even add momentum- 
independent interactions proportional to o\-ro 2 -r to TD to give the more general interaction 

2V = J - (3i(3 2 L + pip 2 C + 751752 if + <J\ • o 2 {-l + fiifoY + p x p 2 T + 751752^) 

+(7i • fa 2 ■ f(N + fcfoT + p lP2 S + 75i75 2 ff), (4.11) 

where 

Pi = A7«, (4.12) 

the covariant a is from Eq. (2.16), f = and N,T, S, and R are arbitrary invariant 

functions of x±. Note that each term in 2V involves identical operators for particle 1 and 
particle 2. As a result, they all commute with each other. For example, we have 

[<7i-r<7 2 -f, 0-i-t7 2 ] = O. (4.13) 

Hence the a single exponential function 2V can also be written as a product of separate 
exponentials. 

Before continuing our discussion on the structure of the covariant Breit equation, we 
consider the remaining linear combination of the two constraint equations (3.3) involving 
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the difference 2(6^ ■ PSi — 9 2 ■ PS 2 ). Using identities in the beginning of this section we find 
that the difference equation becomes 

(ei - e 2 )exp(-V)ij = {H w - H w )exp{V)^ . (4.14) 

Transforming to \& and using the Breit equation (4.10) gives 

, _ (H? Q - HL) T {m\ - ml) T 2(H W + H 20 )P -p T . J , r , 

WW w 

Combined with Eq. (2.13) this gives 

P-p^ = 0, (4.16) 

confirming the expectation that these momenta remain orthogonal after the interaction is 
introduced. This result has also been obtained recently by Mourad and Sazdjian [18] who 
emphasize that this would ensure the Poincaire' invariance of the theory. They further point 
out that the cm. energy dependence of the potential in the "main equation" (our Eq. (4.10)) 
ensures the global charge conjugation symmetry [18] of the system, a feature that is missing 
from the original Breit equation. 

In summary, the constraint equations imply covariant Breit-like equations of a certain 
form (4.10) whose wave function also satisfies the constraint equation (4.16). Alternatively, 
if we start off with a covariant Breit-like equation of the form (4.10) with T> = T>(x±) and 
require simultaneously that P -p^ = 0, we can work backward to obtain the two compatible 
constraint equations (3.3). 

We point out finally that the application of the constraint equations to electromagnetic 
interactions does not involve the term at -fa 2 -f appearing in the Breit interaction Eq. (1.2). 



Nevertheless, it does produce the correct spectrum, as shown in |[il]J12|] using the "external 
potential" form in Eqs. (3.27) of the two-body Dirac equations with L = C = 0, J = —Q. 
As has been recently noted [19], in the context of the Breit-like form Eq. (4.10) of that 
equation, one obtains the reduction 

(#10 + #20 + V\ + V 2ai ■ a 2 + V3751752 + ■a 2 )'V = wV, (4.17) 

which contain pseudovector terms in place of the vector Breit term (a± ■ fa 2 ■ f). 
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V. STRUCTURE OF THE BREIT EQUATION 

The Breit equation (4.10) can be written, as usual, as a set of coupled equations for 
different components of the wave function contained in \1>. We work in the center-of-mass 
system for which P — (1, 0), a — (0, cr), and f = (0, r). We begin by simplifying the general 
interaction (4.11) to the more compact form 

3 

2V = ]T(A, + <n • <t 2 B v + <n r a 2 -rCjqWqW (5.1) 
where the superscripts (1) and (2) label the interacting particles 1 and 2. The operators 

(Qo, Qi, 12, Qs) = (1, 7s, -ip, P) (5.2) 
are defined so that q 1 , q 2 and g 3 are analogous to the Pauli matrices ci, cr 2 and <r 3 satisfying 

QiQj = hi + itijkQk (5.3) 

where i,j and = 1, 2, and 3. Eq. (5.1) is in the form of four-scalar products 

2V = (A + o-i • cr 2 B + a x ■ r <r 2 • r C) ■ Q (5.4) 

involving the "four- vector" Q v = q^q^\ and 

A = (J, H, -C, -L), B = (-/, g, -F, Y), C= (TV, i?, -5, T). (5.5) 

The wave function ^ in Eq. (4.10) can be written as a spinor or column vector with two 
indices, one for each particle 

* = tfWtfC 2 ). (5.6) 

It is however more convenient to express the content of the wave function \1/ in terms of a 
new 4x4 matrix wave function *&' 

where \l/ < ^ 2 - ) has been transposed and an operator a y has been added on the right, as explained 
below. We can represent the operation of for particle 1 and for particle 2 acting 
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on the original spinor wave function \1/ in terms of operations A and B' on the new wave 
function ty', 

AWb^V -> -4^ (1) [^ (2) ] T a^ = AV'B', (5.8) 
where the matrix operator B' is 

S' = a y B T a y . (5.9) 

The arrow in Eq. (5.8) indicates the transformation of the operation of A^B^ on the wave 
function ^ to the operation of A and B' with respect to the new wave function The 
operator a y in Eq. (5.7) insures that operators such as ct^ acting on the second particle is 
represented by 

a (a)# _> ^'^cK T a y = -#'a , (5.10) 

where the same negative sign appears for the different components of the operator a. By 
using the wave function ^' and Eq. (5.8), an operator acting on the first particle appears 
on the left side of the wave function \J r ' ) while an operator O acting on the second particle 
becomes a y O T a y and appears on the right side of 1 J r '. 

In this matrix notation, the righthand side (RHS) of Eq. ( fOIj ) becomes 

(#10 + # 20 )^ ^p-a^' + p-^'a + mi/3*' - m 2 ^' (3 = RHS . (5.11) 

The reduction of Eq. (|4. 10| ) is facilitated by separating the matrix wave function into two 
parts: 

3 

m = ^ + ^ _> = + = £ g K a A ^A, (5.12) 

A,k=0 



where <r = 1, cr = 75 ck 20 



and 



3 

*o = E^o, (5.13a) 

K=0 
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3 3 3 

*i =E^'t = EE Qk^kX- (5. 13b) 

K=0 K=0 A=l 

These parts give different results when operated on by cr x ■ cr 2 , 

OV <r 2 * -> -<r #V = -3*o + (5.14) 

This shows that *o(*o) an d ^i(^'i) are respectively the spin-singlet and spin-triplet parts 
of the wave function. Furthermore, all the operators in Eq. (4.11) commute. With (<Ti • 
r<r 2 ■ r) 2 = 1, the tensor interaction can be written as 

^■vn-TCQ = CQsh ( C .Q) +(Tl -r(T2-r sinh(C • Q) . (5.15) 

This means that interaction (4.11) alone will not mix spin-singlet and spin-triplet states (as 
seen below in Eq.(5.34), the kinetic energy term does mix singlet and triplet states). 

The reduction to the matrix form is easier for the spin-singlet wave function 

= \{\ — (Ti • cr 2 )\I/ because of its simpler spin structure: 

where T = A — 3B — C. A Taylor expansion of the expontial operator on the above equation 
shows that it is necessary to evaluate the basic operation of the type 

T-Q* Q = j^T v qWqW* . (5.17) 

In terms of the new wave function \l/ of Eq. (5.13a), the above equation can be represented 
in the matrix form: 

3 3 3 

H T v (q u q R q' u )^ K Q = ^ T ■ S K q K ^ K0 , (5.18) 

where q' u = a y q^a y = e u q u , with (e , ei, e 2 , e 3 ) = (1,1,1,-1) and we have introduce the 
quantity (S K ) V defined by 

QvQkcL = (S K )uq K , (5.19a) 

and 
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{3 K ) V = e„[l + 2(1 - 6 K 0)(1 - 5 u0 )(S KU - 1)] . (5.19b) 

The 4x4 matrix S is called a signature matrix because its matrix elements can only be +1 
or — 1 . The row vectors of S are 

5 = (1, 1, 1,-1), 

51 = (l, 1,-1, 1), 

5 2 = (l,-1, 1, 1), 

5 3 = (1,-1, -1,-1). (5.20) 
Eq. (5.18) can be applied repeatedly to give the desired result 

3 

-> w e^q^o = LHS . (5.21) 

It can also be used to prove the general result 

f(T ■ Q)^ - X! f(T ■ S K )q K ^ K0 . (5.22) 

The treatment of the spin-triplet expression 

e 2V m 1 = e {A+B) Q [cosh(C ■ Q) + a l ■ r <r 2 • r sinh(C ■ Q)]^> 1 (5.23) 

is simplified by noting that the q K and er matrices are independent of each other. Hence the 
Q dependences can be eliminated in favor of the signature vector S K with the help of Eq. 
( |5.19a| ), which also applies to the q structure of This leaves the spin-dependent part 
which has the form 

(j\ ■ r cr 2 ■ r^i — > — cr ■ f^/cr ■ r 

= Y j q K {<r-^ K -2 ( T-rr-^ K ). (5.24) 

Hence, using Eq. ( |5.15| ), the spin-triplet part of the left hand side of Eq. (4.10) is 
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LHSi = w e {A+B) - SK q K [e c - SK (T -i/) K -2 sinh(C • S k )(t -rr-^ K }. (5.25) 

K=0 

The Breit equation (4.10) for the matrix wave function is thus 

LHS + LHSi = RHS (5.26) 

where the expressions are from Eqs. (5.11), (5.21), and (5.25). 
Eq. (5.26) can be written explicitly as 

3 f 3 1 

Y Qk(T\\ W kX ^k\ - (1 - Sxo)V K r x Y r^nv \ 

3 r 3 r i 

= Y Ys\ a i^ a >~ + 9«^AQ!i \Vi + m^q^x - m 2 q K a x q3 

K X=0 L i=l ' 

where 



Ax , (5.27) 



1^ kA = w e ( A - 3B - c > s «5 X0 +w(l- 5 X0 )e^ A+B+c ^ , (5.28a) 

% = 2we {A+B) - s * sinh(C • S K ) . (5.28b) 

Multiplying the equation from the right by <r x q K and taking traces we finally get a set of 16 
coupled equations for the wave function components, 

\ W KX ip K \5\ - (1 - Sxo)V K f\ Y \ 

3 3 3 

= Y Yi 1 + fK'ifx'i)g K iK'9xix'Pi^K'X' + Y ( m i - f K '3m 2 )g K 3 K '^ K 'x , (5.29) 

k'X'=0 i=l k'=0 

where 

f Kl = 5 K0 + (l-5 K0 )(25 Kt -l), (5.30) 

and 

9kik' = 5 K0 5 iK/ + 5 K ' 5 iK + (1 - S k0 )(1 - 5 K/0 )ie KiK , . (5.31) 
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The structure of these equations becomes more transparent by writing them out explicity 
in terms of the singlet and triplet wave functions 

Ooo, ^10, ^20, ^so) = O, 0, ix, V), (5-32) 
(^o, -01, "02, ^s) = (*!>, <t>, iX, V) ■ (5-33) 



where 



[w - U^ifj = 2p ■ 4> + (mi - m 2 )r), (5.34a) 

[w - U^](f) = 2p ■ %l> + (mi + m 2 ) X , (5.34b) 

[w - U x ]x = (mi + m 2 )(f), (5.34c) 

[w - U v ]rj = (mi - m 2 )i> ) (5.34d) 

[w — V^- V^rr-]ij) = 2p<f) + (mi - m 2 )rj, (5.34e) 

[w-Vtj,- V<prr-]<f> = 2pip + (mi + m 2 )x, (5.34f) 

[w — V x — V x ff-]x = —2ip x 77 + (mi + m 2 )(f>, (5.34g) 

[w — V v — V v rr-]r) = —2ip x x + ( m i — m 2 )tp, (5.34h) 



w 



w - V K = we( A+B+c > s «, (5.35b) 
and V K is defined in Eq. (5.28). 
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In order to see explicitly the distinction between the traditional Breit approach and that 
of constraint dynamics, we perform an angular momentum decomposition. For the spin- zero 
wave functions we take 



i> = ^jYjmi = (pjY jm , x = XjYjm, V = TjjY jm , (5.36a) 

where Yj m is an ordinary spherical harmonic (j = I here). For the spin-one wave functions, 
we take a form that depends on the spatial parity: 

= <t>jy Y j™ + <f>jzZ jm , or = (f) jx X jm ; (5.37a) 



if) = iJjyY jm + 4> jz Z jrn , or = ip jx X jrn ] (5.37b) 

^ = VjyYjm + VjzZjm, or = Vjx^jm'j (5.37c) 

X = XjyYjm + XjzZj m , or = XjxXj m ] (5.37d) 

where the three orthonormal vector spherical harmonics are 

L r rp 

Xj m / = Yj m , Yj m Yj m , Z jm i i = Yj m . (5.38) 

The first and last vanish for j = 0. 

To obtain the radial wave equations, we use the following identities 

v 

- ■ X jm = 0, p ■ X jm = 0, (5.39a) 

r Y .„ v >Jj(3 + l) v Z jm 

— X Ji-j m — lZij m , p X A-j m — - Ijm H ~ , (O.oyDJ 

r 1% 

~ ' Yj m Yj m , p ■ Yj m —Yj m) (5.39c) 



r Jj(j + 1) 

- x Y jm = 0, p x Y jm = -X jm , 5.39d 

r r 
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r + 1) 

- • Z jm = 0, p ■ Z jm = i± Y jm , (5.39e) 

r r 

- x Z jm = X jm , p x Z jm = (5.39f) 

r r 

The equations so obtained can be divided into two sets of different total parity. One set has 
the natural parity solution: 

V* jm = (-) j+1 y jm , (5.40) 
and involves the 8 wave functions: 

= (/>jY jm , x = XjYjm, (5.41a) 

4> = (f) jx X jm , tj) = iip jy Y jrn + iij jz Z jm , (5.41b) 

"X. Xjx-X-j m , tj iT]jyYj m + if]j Z Zj m . (5.41c) 

(Four of the wave functions ^jxi^jziXjx an d m 2 do not appear for j = 0, because their 
angular parts vanish.) The eight simultaneous equations satisfied by them are 

4 Jj(j + 1) 

(w - U<t)<j>j = 2ip' jy + -i/j jy - 2-i — 4> jz + (mi + m 2 )xj, (5.42a) 

(w -V^- V^ipjy = -2<j)j + (m 1 - m 2 )i] jy , (5.42b) 

2 

(w - V v )r] jz = -2x' jx - -Xjx + (mi - m 2 )^ jZ) (5.42c) 

(w - V x )xj X = 2rj' jz + %ij g - 2 ^^^ \ jy + (mi + m 2 )(j) jx . (5.42d) 

(w - U x )xj = (mi + m 2 )^, (5.43a) 



( «' - V^,)-0j Z = -2- 0j + (mi - m 2 )r) jz , (5.43b) 
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(w - V^)(j> jx = (mi + m 2 )xjx, (5.43c) 
Jj(j + 1) 

(w-V v - V v )r] jy = -2 Xjx + (mi - m 2 )4) jy . (5.43d) 

The second set has the "unnatural" parity 

V* jm = (-) j * jm , (5.44) 
and involves the 8 wave functions: 

ip = iJjY jm , 7] = mY jm , (5.45a) 

-0 = i/jj X X jm , 4> = i(f) jy Y jm + i(f) jz Z jm , (5.45b) 

"H Vjx-^-jm, X ^XjyYjm "i" ^Xjz^jm- (5.45c) 

The eight simultaneous equations satisfied by them are 

4 Jj(j + 1) 

(w - U^)i)j = 2<p jy + -<p jy - 2- <p jz + (mi - m 2 )m, (5.46a) 

2 

(w - ^x)Xj 2 = — 2?7j- x - -m x + (mi + m 2 )0 i2 , (5.46b) 

(iw - K,)^ = 2Xj- 2 + ^Xj 2 - 2 ^ 3 ^ r + ^ Xjy + (™i - rn 2 )^ jx . (5.46c) 

(iu - - V^jy = -2^ + (mi + m 2 )xjy, (5.46d) 

(w - Ur,)rfj = (mi - m 2 )ipj, (5.47a) 

(w - V^)ip jx = (mi - m 2 )r} jx , (5.47b) 
Jj(j + 1) 

(tw - V^,)0^ = -2-i + ( m i + m 2 )xj Z , (5.47c) 
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Jj(j + 1) 

(w-V x - V x )xjy = rj jx + (mi + m 2 )<p jy , (5.47d) 

where the terms w — U, w — V are given in Eqs. (5.35a-b) and w — V — V by 

w-V K -V K = we (A+B - c) - SK . (5.48) 

Each set of 8 equations consists of 4 differential equations [Eqs. (5.42) or (5.46)] and 4 
algebraic equations [Eqs. (5.43) or (5.47)]. The algebraic equations can be used to express 
4 of the 8 wave functions in terms of the remaining 4. For example, two of the eliminated 
wave functions are 

* = ^*, 

* = ^. (5.50, 

After the elimination, pole singularities could appear in the differential equations at particle 
separations r where equations such as Eq. (5.49) have poles, e.g., where 

w — U x = 0; w — U v — 0, mi ^ m 2 (5.51) 

provided that the total relativistic center-of-mass energy w is nonzero. These are the well- 
known singularities that plague the traditional Breit equations. 

However, from the perspective of the exponentiated interactions of constraint dynamics 
such as that shown in Eq. (5.35), these structural poles can appear only if the potential 
generators in the exponent go to — oo. In the absence of such singular behavior, the constraint 
two-body equations, or their Breit analogs, are free of the pathologies described in [6,7]. 
These pathologies arise because the wave functions appearing on the right-hand side of Eq. 
(5.49) must have zeros at the pole positions. These additional boundary conditions give rise 
to spurious resonances in the continuum for sufficiently strong interactions, and to spurious 
bound states for any nonzero interaction strength at total energies (including rest masses) 
which go to zero. 
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We should add, for the sake of completeness, that another class of spurious solutions 
appears when the values of the exponential generators go to +00. Then the algebraic 
equations impose the new boundary conditions that the wave functions appearing on the 
left-hand sides of these equations must vanish at these singular points. 

It is important to point out that if the exponentiated potentials given in Eqs. (5.35) and 
(5.48) are approximated by finite Taylor approximants, it might be possible to satisfy the 
singularity conditions such as (5.50) whenever the approximant has a zero where the true 
value of the exponential function is nonzero. In other words, it is the exponential structure 
of the effective interaction that protects the CTB Dirac equations from the undesirable 
pathology. We trace that exponential structure to the hyperbolic forms (3.3) of the constraint 
equations which in turn are motivated by the requirements of compatibility. Recall that in 
those forms the "generators" A produce potential terms in the "external potential" form of 
the constraint equations which have the expected transformation properties, including the 
essential recoil corrections. 

Our result can also be interpreted in another way. The exponential functions appearing 
on the left hand side of Eq. (4.10) can be expanded in terms of hyperbolic functions. Indeed, 
one could express the "effective potentials" contained in exp(2V) in the form given by the 
right-hand side of Eq. (4.11). If we now parametrize these effective potentials directly using 
functions with no singularity at finite r, we will find that under favorable circumstances the 
singularity condition (5.50) can still be satisfied. This is in fact the result of [6,7]. Thus the 
regularity of the effective potentials appearing in exp(2T>) does not guarantee the regularity 
of the resulting Breit equation. It is the regularity of the basic potentials appearing in the 
exponent 2V that guarantees the regularity of the resulting Breit equation. 
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VI. COMPARISON OF THE CONSTRAINT AND BREIT EQUATIONS FOR 

QED 



In this section we discuss the implications of our new Breit-like equation in quantum 
electrodynamics. Consider first the original Breit equation whose matrix form is 

(w + — W + a- a ■ ^>'ct + b—OL ■ r^'ot f = p- a*' + p ■ ^'a. + - m 2 ^'(3 , (6.1) 

r r r 

where a and b are potential parameters. By identifying it with the center-of-mass form 
of the Breit-like form Eq. (4.10) of our covariant constraint equation, we can solve for the 
generators given by Eq. (5.5). Using Eqs. (5.26) or equivalently Eqs. (5.27-28) we obtain 
the following twelve equations for the potentials. 




(6.2a) 




(6.2b) 



exp[(A 



3B 



C)-S 2 ] = l + (1 



3a - 6)C, 



(6.2c) 




(6.2d) 



exp[(A + B + C)-S ] = 1 + (1 



a - b)(, 



(6.2e) 



exp[(A + B + C) ■ S x ] = 1 + (1 



a - 6)C, 



(6.2f) 




(6.2g) 



exp[(A + B + C) ■ S3] = 1 + (1 + a + 6)C, 



(6.2h) 




(6.2i) 
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exp[(A + B-C)-S 1 ] = l + (l-a + b)(, 



exp[(A + B-C)-S 2 ] = l + (l + a- b)(, 



(6.2k) 



exp[(A + B - C) ■ S 3 ] = 1 + (1 + a - b)(, 



(6.21) 



in which ( = a/(wr). The first four equations come from the singlet wave function (A = 
terms in Eq. (5,28a)), while the last eight arise from the triplet, with the last four coming 
from a combination of Eq. (5.28b) and the triplet (or A ^ 0) part of Eq. (5.28a). These 
twelve algebraic equations can be solved for the twelve unknown generators shown in Eq. 
(5.5). One can readily show that six of these generators vanish (corresponding to scalar, 
pseudoscalar, and tensor interactions): 



The remaining generators, J, H, J, G, iV and R, are nonzero, but we shall not need them in 
our discussion except in the combinations appearing on the left-hand sides of the algebraic 
equations, Eqs. (5.43) and (5.47). These are just the six equations shown in Eqs. (6.2c-f) 
and (6.2k,l), now expressible as 



w-U x = w-U T1 = w exp[J + 31 — N — H + 3Q + R] = 1 + (1 — 3a — &)( , (6.4) 



w -V x -V x = w -V v -V r) = w exp[J -I + N-H + g-R] = l + (l + a-b)(. (6.6) 



Singularities at finite particle separation r arise, for nonzero w, when the right-hand sides 
of these equations vanish (assuming that the numerators in the algebraic equations in which 
these occur do not). Recalling the context in which these equations appear, we see that this 
occurs when 



S = T = T = Y = C = L = 0. 



(6.3) 



w-V i , = w-V 4> = w exp[J -I + N + H + g + R] = l + (l-a-b)(, 



(6.5) 



3a + b > 1 



(6.7) 
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for all V = (— states, and all V = states when m 1 ^ m 2 ; 

a + b>l (6.8) 

for all V = states, all V = (— states when mi 7^ m 2 ; and for j = states when 
m i = m 2 an d 

b-a>l (6.9) 

for all V = (— Y +1 states when m x 7^ m 2 , and for j = states when mi = m 2 and all 
V = (~y states. This confirms the result first derived in Ref. [7]. 

It is worth noting that the original Breit equation corresponds to a — b — 1/2 while 
the Barut equation corresponds to a = 1,6 = [21]. Both equations are therefore singular 
according to the first of the above three conditions. 

From the perspective of the exponential generators of constraint dynamics, the vanishing 
of the right-hand sides of Eqs. (6.4-6.6) is possible only when the generators take on the 
value —00 at finite r. Hence they can be avoided by simply using generators not having 
such r singularities. Indeed, in their constraint approach to QED, Crater and Van Alstine 
[11,12] found that 

e G = e' J = (1 + 2a/{wr)y 1/2 (6.10) 

with all of the remaining ten generators set equal to zero. As a result, no midrange zero 
appears since 

w — U Y = w exp(J + 3£) = w exp(2£) = ; — (6.11) 

x l + 2a/wr 

is finite for positive r. Furthermore, the Lorentz nature of the two nonzero generators 
matches that of the vector character of the QED interaction at this level. In contrast, the 
Breit generators involve pseudovector as well as vector parts and additional vector and pseu- 
dovector "tensor terms" corresponding to the Coulomb gauge ( as opposed to the Feynman 
gauge used in the constraint equations). 
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The conclusions obtained in this section for QED can be extended to other interactions 
(scalar, pseudoscalar, psuedovector and tensor) important in semiphenomenological appli- 
cations in nuclear and particle physics. Thus, one may with safety solve the new Breit 
equations numerically with no concern about structural singularities which would otherwise 
render such solutions meaningless. 

VII. CONCLUSIONS AND DISCUSSIONS 

Breit equations and constraint Dirac equations of relativistic two-body quantum mechan- 
ics have markedly different properties: Breit equations could have structural singularities at 
finite particle separations even when the interactions themselves are nonsingular there, while 
constraint Dirac equations seem to be free of them. CTB Dirac equations are manifestly 
covariant, while Breit equations are not, being valid only in the center-of-mass frame. We 
are able to understand, and to reconcile, their differences in this paper. 

The constraint Dirac equations were originally derived for scalar interactions with the 
help of supersymmetries in addition to Dirac's constraint dynamics. Generalizing this con- 
cept to arbitrary interactions leads to a "hyperbolic" form of the equations. The two 
single-particle Dirac equations can then be recast into two other equations: (1) a covariant 
Breit-like equation with exponentiated interactions, and (2) a covariant equation describing 
an additional orthogonality condition on the relative momentum. We use the equivalence 
between these two types of equations to show that the constraint Dirac equations are com- 
pletely free of the unphysical structural singularities when the exponential structure of the 
interactions are not tampered with. 

The advantage of the constraint form of the Breit equation is that the structural singular- 
ities of the traditional Breit equation are now moved entirely to the exponential generators 
of the interaction. As a consequence, they can be eliminated right from the beginning by 
the simple requirement that these generators themselves be free of singularities at finite 
separations. The resulting Breit equations are then guaranteed to be free of the undesir- 
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able structural singularities that plague traditional Breit equations. These improved Breit 
equations, which are dynamically equivalent to the constraint Dirac equations, can now be 
used in nonperturbative descriptions of highly relativistic and strong-field problems such 
as those appearing in two-body problems in quantum electro- and chromo-dynamics, an in 
nucleon-nucleon scattering. 

Of course, the constraint Dirac equations can also be used, now that we know how to 
keep them singularity-free. However, in actual applications, they have to be reduced down 
to a set of coupled differential equations before actual solutions can be attempted. These 
differential equations are transforms (see Eq.(4.9)) of those obtained from the equivalent 
Breit equations. So in reality the two different formulations have now become completely 
identical to each other. 
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